Galileo satellite constellation and extensions to General Relativity 
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We consider the impact of some known extensions of General Relativity in observables that will 
be available with the Galileo positioning systems, and draw conclusions as to the possibility of 
measuring them. We specifically address the effects of the presence of a cosmological constant, a 
Yukawa-like addition to the Newtonian potential, and the existence of an extra, constant accelera- 
tion. We also consider the phenomenological impact of a broad class of metric theories, which can 
be classified through the parameterised Post-Newtonian formalism. 
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I. INTRODUCTION 



Q ' The Galileo positioning system poses a great opportunity, not only for the improvement and development of 
, new applications in navigation monitoring and related topics, but also possibly for fundamental research in physics. 
Indeed, together with the already deployed american and russian counterparts, the Global Positioning System (GPS) 
and Glonass, satellite navigation may be considered the first practical application where relativistic effects are taken 
into account, not from an experimental point of view, but as a regular engineering constraint on the overall design 
requirements. Indeed, effects arising from special and General Relativity (GR) - gravitational blueshift, time dilation 
^ and Sagnac effect - may account to as much as ~ 40 ^s/day, which is many orders of magnitude above the accuracy of 
Y the onboard clock deployed in these systems. Moreover, the gravitational Doppler effect, of the order of VV/c^ ~ 10~^° 
^ . (where Vjv = GMe/Re is the Newtonian potential, G is Newton's constant. Me ~ 6.0 x 10^** kg is the Earth's mass. 
Re ~ 6.4 X 10^ m is its radius and c is the speed of light) falls within the 10~^^ frequency accuracy of current 
space-certified clocks, and must therefore be taken into account: in GPS, this is done by imposing an offset in the 
i-H ■ onboard clock frequency, while in Galileo this correction should be corrected by the receiver. For further details, the 
" reader is directed to Refs. [1-4] and references within. 

This said, it is not clear as to what extent the accuracy of the Galileo positioning system may be improved - which 
is designed to offer pinpoint localisation within an error margin of 1 to, against the 10 to margin of previous the 
, GPS system - so to provide clues to the nature of models beyond the current GR scenario. In this study, we aim at 
• ' establishing some bounds on the detectability of commonly considered extensions to GR [5]. This paper is organised 
^—^ , as follows: firstly, we assess the main relativistic effects that are present in the Galileo system. We proceed and 
consider the possibility of measuring a variety of extensions of GR and conclusions are then drawn. 



II. MAIN RELATIVISTIC EFFECTS 



$_( ' A. Frame of reference 

Assuming that time-dependent effects are of cosmological origin, and hence of order where Hq is Hubble's 

constant, one may discard these as too small within the timeframe of interest; hence, one assumes a static, spherically 
symmetric scenario, posited by the standard Scharzschild metric. In isotropic form, this is given by the line element 



ds^ = - (l + ^) (c dtf + Y^^dV + (c dtf + - ^) dV , (1) 

where dV — dr^ + dft^ is the volume element, and V is the gravitational potential. In the standard GR scenario, the 
latter coincides with the Newtonian potential V = Vn — ~GME/r{\ -I- E"^]^ J„), where the J„ multipoles account for 
the effect of geographic perturbations and density profiles. 
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However, one must introduce the rotation of the Earth with respect to this fixed-axis reference frame, with angular 



velocity w = 7.29 x 10 ^ rad/s; by doing a coordinate shift t' 
Langevin metric, given by the line element 



t r' 



e' = 9 and 



Lot', one gets the 



2V ( u)rsmi 
V c 



(c dtf + 2wr^ sin^ 6'#di + ( 1 + ^ ) 



(2) 



where, for simplicity, primes were dropped. Asides from a non-diagonal element, one obtains an addition to the 
gravitational potential, which could be viewed as a centrifugal contribution due to the rotation of the reference frame. 
One can then define an effective potential $ = 2V — (uirsinO)'^; the parameterization of the Earth's geoid is obtained 
by taking the multipole expansion of V up to the desired order and finding the equipotential lines $ = $o (the latter 
being the value of $ at the equator), and solving for r{0, </>). 

In the above line elements, the coordinate time coincides with the proper time of an observer at infinity. However, 
since one wishes to evaluate the ground to orbit clock synchronisation, it is advantageous to rewrite the metric in 
terms of a rcscaled time coordinate, which coincides with the proper time of clocks at rest on the surface of the Earth; 
this is best implemented by resorting to the above-mentioned geoid, since its definition as an equipotential surface 
$ = $0 indicates that all clocks at rest in it beat at the same rate; hence, rescaling the time coordinate according to 
t ^ {1 + $o/c^)ii one gets the metric given by the line element 



ds' 



2($-$o) 



2V 

(c dty + 2ur'-sin^9d<pdt [l - — ]dn . 



(3) 



Finally, if one reassumes a non-rotating frame, the metric is given by the line element 



ds^ 



1 + 



2(y-$o) 



(c dtf + 1 



2y 



dn 



(4) 



B. Constant and periodic clock deviation 

One may now consider the difference between the time elapsed on the ground and the satellite clock; keeping only 
terms of order c~^, one finds that the proper time increment on the moving clock is approximately given by 



V -^0 v'^ 
dr = ds/c = ( H ^ 1 dt . 



(5) 



Considering an elliptic orbit with semi-major axis a, and taking V = Vn » GMs/r, this may be recast into the form 

[1] 



dr = ds/c 



3GMe <^>o 2GMe 
2a(? (? (? 



ii- 









dt 



(6) 



The first constant rate correction terms in the above amount to 



3GMe <E>o 



2ac2 



= -4.7454 X 10 



-10 



(7) 



for the Galileo system, and —4.4647 x 10""'^'', for the GPS system; this indicates that the orbiting clock is beating 
faster, by about 41 jis/day, for the Galileo system, and 39 iis/day, for the GPS system . For this reason, the GPS 
system has a built in frequency offset of this magnitude, while the increased computational capabilities made available 
to current and future receivers of the Galileo system leave this correction to the user. The residual periodic corrections, 
proportional to 1/r — 1/a, have an amplitude of order 49 ns/day, for the Galileo system, and 46 ns/day, for the GPS 
system. 
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C. Shapiro time delay and the Sagnac effect 

The so-called Shapiro time delay, a second order relativistic effect due to the signal propagation is given by [1] 

^ol ^GMe , a I \ 

Ataelay = ^ + ^^ln[l + —j , (8) 

where we have integrated over a straight line path of (proper) length I. Evaluating this delay, one concludes that this 
effect amounts to 6.67 x 10~^^ s. 

Also, one must consider the so-called Sagnac effect, which arises from the difference between the gravitational 
potential V and the effective potential when proceeding from a non-rotational to a rotational frame. Hence, one 
gets tlic additional time delay 



Atsagnac = ^ / d(j) = — dA^ , (9) 

J path ^ J path 

where dAz is the orto-cquatorial projection of the area clement swept by a vector from the rotation axis to the satellite. 
For the Galileo system, this yields a maximum value of 153 ns while, for the GPS system, one gets 133 ns. 

One concludes this section by recalling the main effects affecting the considered global positioning systems: a 
frequency shift of order 10^"'^*' and a propagation time delay (Shapiro plus Sagnac effect) of the order 10~^ s. In what 
follows, one shall compute the additional frequency shift and propagation time delay induced by common proposals 
for extensions of GR, and compare the results with the above quantities, plus the frequency accuracy of 10~^^ and 
the time accuracy of Galileo, of order 10~^ s, which corresponds to a optimistic spatial accuracy of 30 cm. 



D. Post-Newtonian effects 



We address now the issue of measuring Post-Newtonian effects with the Galileo positioning system. As the moniker 
indicates, these are effects below the Newtonian order, that is, GMe/Rec^ ~ 10~^°. A general formalism exists 
to describe lower-order effects induced by extensions to GR and alternate theories that resort to a metric approach 
of gravity; indeed, any such theory may be analysed locally and compared with the so-called Parameterised Post- 
Newtonian (PPN) metric [6, 7], given by the line element 



2V (V 



(c dtf + 



V 

1-27^ 



dV 



(10) 



In the above, the parameter j3 measures the non-linearity of the superposition law for gravity, while 7 indicates the 
space curvature produced per unit mass. For clarity, we consider only a simplified version of the full PPN metric; the 
latter encompasses ten PPN parameters, characterising the underlying fundamental theory, and may be related to 
violation of momentum conservation, existence of a privileged reference frame, amongst others. GR is characterised 
by P = 7 = 1, while all remaining parameters vanish. For that reason, the quantities P — 1 and 7 — 1 measure the 
deviation from the predictions of the currently accepted theory. Experimentally, it is found that — 1| < 2 — 3x 10"'* 
(Nordtvedt effect) and 7 - 1 = (2.1 ± 2.3) x 10~^ (Cassini radiometry). 

Unfortunately, it is clear from the above equation that Post-Newtonian effects arise only at an order ~ 10" 
undetectable by the accuracy of the GPS and Galileo systems. 
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III. DETECTION OF THE COSMOLOGICAL CONSTANT 



Latest observations indicate that the Universe is experiencing an accelerated expansion, which may be characterised 
by the presence of a cosmological constant A ^ 10^''^ s^^, acting as a negative-pressure fluid (see e.g. [8] and 
references therein). By matching the outer Friedmann- Robertson- Walker metric with a static, symmetric solution 
given by Birkhoff 's theorem, we may derive the Schwarschild-de Sitter metric, given by the line element (in anisotropic 

form) [9], 
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This indicates that the cosmological constant induces an additional term to the potential, of the form Va — — Ar^/6; 
since its expected effect is assumed to be small, one may neglect the issue of performing a coordinate change to an 
isotropic, co-rotating frame of reference, as well as the identification of proper time with clocks at rest on the surface 
of the gcoid (however, notice that the identification of proper time as that measured by a clock at rest at infinity 
breaks down, due to the Schwarschild "bubble" breaking down at a distance r given by the condition V„ = Va). 

The frequency shift of a signal emitted at a distance from the origin r = Re + h (for the Galileo system, h = 
17.2 X 10^ m) and received at a distance r = Re is given by 



fEarth\ _ jgm Earth _ / \-2V{Re)Ic'^ V {Re) ' V {Re + h) 

12 • v-"-^; 



fSat J V 500 Sat 2V{Re + h) / C2 

Hence, one may compute the additional frequency shift induced by this extra potential contribution, through 

fEarth\ ^, V^{Re) - V^Re + h) _ A ^^^i?^ ^ ^ ^^.38 ^ (^3) 



fsat /a c2 6c2 

which clearly falls bellow the accuracy tf^ = 10"^^ of the Galileo constellation. 
Likewise, the propagational time delay is given by 



j-RE+h 
' Re 



1 rrtE-ta 

- / V(r) dr . (14) 

cJre 

Hence, the cosmological constant induces a further delay of 



1 pRB+h » 2 K 

^*^^cj dr=j^h[{3RE{RE + h) + h^]^ 10-^° s , (15) 

also many orders of magnitude below the time resolution of 10^^ s. Therefore, one concludes that the cosmological 
constant is completely undetectable by the Galileo system. 



IV. DETECTION OF ANOMALOUS, CONSTANT ACCELERATION 

An anomalous constant acceleration could model first-order effects arising from some fundamental theory of gravita- 
tion which expands upon GR, or indicate some threshold between known dynamics and yet undetected, exotic physics. 
One examples stems from the so-called Modified Newtonian Dynamics (MOND) model [10-12], which attempts to 
account for the missing matter problem indicated by galactic rotation curves without the need for dark matter, 
by featuring a departure from the classical Poisson equation at a characteristic acceleration scale of 10^^*^ mj . 
Also, although yet unmodelled or theoretically unaccounted for, an anomalous, sunbound, constant acceleration 
a = (8.74 ± 1.33) x 10"^° m/s^ has been reported to affect the Pioneer 10/11 probes [13-15]. 

An anomalous, constant acceleration a may be phenomenologically modelled by a potential Va = ar; following the 
procedure depicted in the previous section, the following frequency shift is obtained 

(fEarth\ Va{RE)-Va{RE + h) ^ oh 

V fSat C2 C2 • 

Comparing with the frequency accuracy e/^ = 10"^^, one finds that only accelerations a > (?ef^/h 10 
be detected. 

The propagational time delay due to this extra potential addition is given by 
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1 r^'^'^+'^ ar a 
^ta = - — dr= ^h{2RE + h) , (17) 

C J Re C IC 

and comparison with a time accuracy of 10~^ s yields the condition for detectabihty a > 100 m/s^. Therefore, one 
concludes that accelerations of the order 10~^° — 10~^ m/s^ are beyond the observable reach of the Galileo system; 

conversely, dctcctability of a constant acceleration of the order of 10"^*^ m/s^ would require an improvement of 7 
orders of magnitude in frequency accuracy (to e/^ ~ 10^^''') and 12 orders of magnitude in time resolution (to 10"'^^ s). 



V. DETECTION OF YUKAWA POTENTIAL 



A common phenomenological approach to extensions of GR lies in assuming that the potential has a modified 

Yukawa form, 



y(r)=-i^(l + ae-^A) , (18) 

where a is the strength of the perturbation, A its characteristic range, and Goo the gravitational coupling for r ^ oo; 
the latter may be regarded as a redefinition of Newton's constant G, through G = Goo(l + a)- This potential may 
be separated into a Newtonian-like potential and an extra potential Vy = —{aGME/{l + a)r)e~^^^\ The Yukawa 
contribution may arise from scalar/tensor field models, where the range is related to the mass m of the scalar field, 
A oc [5]. 

Tight experimental constraintK are available, stemming from sc;vcral sources and regimes, as may be seen in Fig. 
1. Clearly, two yet unexplored avenues remain: the sub-millimeter regime, A < 10""^ m [16], and an astronomical 
regime, A > 10^^ m » 0.1 ly. 

Following the previous steps, one first obtains the extra frequency shift 



f Earth \ 

^ Vy{Re) - Vy{Re + h) _ GMe f a \ ( u,^ Re 



fSat Jy ^ C2i?£; Vl + ay V ^ 

The additional time delay is given by 
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The above expressions may be considerably shortened if it is assumed that this additional "fifth-force" is a long-range, 

A ^ r, or short-range interaction, A <C r. 

A. Short-range fifth force 

If the range of the Yukawa interaction is short-ranged, A <C /i. Re, one obtains 

' fEarth\ GMe f a \ ^-Re/\ _ (21) 



fsat Jy c^iiis \l + a 
If this effect is undetectable within the frequency accuracy e/^, one obtains the constraint for small a 



GMe 
c^Re 



gfle/Ag^^ w 1.4 X IQ-^e^"/^ > 1 , (22) 



which yields no new insight into the yet unexplored sub-millimetric regime, as depicted in Fig. 1. 
Likewise, the additional propagation time delay is given by 

My = ^ + — \ , (23) 



(? \ Re 

so that comparison with the time accuracy of At = 10~^ yields, for a <C 1 



a < 



At« 50 . (24) 



Hence, one concludes that the short-range regime of a hypothetical Yukawa fifth force cannot be probed by the 
Galileo system. 

B. Long-range fifth force 

If one follows the inverse assumption of the previous subsection, and assumes a long range fifth force, A 3> h,RE, 
the exponential terms may be expanded to first order in r/A; as a result, the induced propagation time delay becomes 

GMeu h 

^tY^ ■ (25) 

If the effect is undetected at a level of accuracy At ~ 10~^ s, one obtains, for small a 

|„|<f!^^«4xl0-«f-^') . (26) 
' ' GMe h \lmj ^ ^ 

For a lower bound of A w 10^ m (only one order of magnitude above RE,h), we obtain the incompatible result 
a < 400. 

Regarding the additional frequency shift, one finds 



f Earth \ ^ GMsa h 

fSat )y~ 2c2A' 



^ - (27) 



so that comparison with the frequency accuracy level of e/^ ~ 10 yields, for a <C 1 



or, equivalently, a quite interesting bound 



log a < -5 + log ef^ +2 log ^-^^^ . (29) 

One may plot the different constraints obtained by varying the frequency accuracy e/^, as seen in Fig. 1; this shows 
that, at the current level, no new bounds are produced. Also, it shows that, at a level e/^ ~ 10~^^, the region below 
the "trough" at A ~ 10^ m (corresponding to a < 10~^) could be investigated. 



VI. CONCLUSIONS 



In this work, we have addressed the possibility of detecting signals of new physics through the use of the Galileo 
positioning system. This application could be valuable, as any unexpected new phenomenology could provide further 
insight into what lies beyond General Relativity. We have specifically looked at the propagation time delay and 
frequency shift induced by three different models, namely a potential related to the presence of the cosmological 
constant, the influence of an anomalous, constant acceleration, and the addition of a Yukawa- like fifth force. We also 
briefly discussed the (im) possibility of measuring post-Newtonian effects with the Galileo system. 

Unfortunately, our conclusions indicate that the available observables are not suitable for the intended purpose; 
indeed, while these render the detection of the cosmological constant unattainable, they also indicate that the current 
accuracy is many orders of magnitude above that needed to probe interesting regimes of anomalous constant accel- 
eration (a ~ 10~^° — 10~^ m/s^) or Yukawa range A > 10* m and strength a < 10^*. Indeed, a frequency accuracy 
of lO"'^^, near the "quantum" regime, is required to further probe the desired scales. Although this is clearly beyond 
the grasp of any foreseeable global positioning systems, and yet unavailable in space certified clocks, such precision 
might be attainable in the future. 

Finally, we remark that, although it was not the purpose of this study, the Galileo positioning system could be 
paramount in improving the bound on violation of the Local Positioning Invariance (LPI) principle [5]; this tenant, 
one of the fundamental pillars of General Relativity, postulates that clock rates arc independent of their spacetime 
positions. Experimental constraints on allowed relative frequency deviations indicate that this invariance holds down 
to a level of 2.1 x 10~^ [17]. Endowing one or more elements of the Galileo constellation with higher precision clocks 
and allowing for sufficiently stable communication with stations on Earth, possibly through a microwave link, could 
yield an improvement of up to two orders of magnitude on the LPI. Another alternative could involve installing 
cornercubes on the surface of one or more elements of the Galileo system, so to allow for accurate laser ranging. It is 
tempting to call this subset of the Galileo constellation Siderius Nuncius, the Celestial Messenger, given its potential 
in helping to sort out the mysteries of the Cosmos. 
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